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ABSTRACT 


The problem of analyzing time limited complex wave forms 
having time variant frequency domain characteristics is discussed. 
A bell tone is selected as a wave form to analyze and it is then 
synthesized to produce an approximation to the original sound. 


An electronic device is constructed to simulate all required fog 


signals for a sailboat, including a rapidly ringing bell. 
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SECTION 1 


INTRODUCTION 


The primary objective of this thesis has been to harmonically 
analyze a complex wave form, and then.to synthesize this tone using 
solid state circuitry. The sound of a.ringing bell was chosen to 
be evaluated since it represented the most.difficult class of wave 
forms to analyze. The main form of analysis was repetitive sliding 
time windows of discrete data which were transformed.to.the frequency 


domain. A Fast Fourier Transform algorithm was used to transform the 


.discrete data. ‘The techniques are not original with the author, but 


they represent a relatively new application of discrete Fourier 
analysis on a general purpose digital.computer. This method.of 
analysis is applicable to any discipline wherein frequency spectrum 
information is desired. Recent and future projects at.the U..S. Naval 
Postgraduate School in this.area include, but are not limited.to, 


voice pattern recognition, helium atmosphere voice distortion, the 


.study of surface waves on water, and squirrel heart-rates under 


stimuli. 

As a second method of spectral analysis, an analog narrow band 
spectrum.analyzer was employed-to check the results of .the discrete 
method. The.comparison.of results was .favorable.and.is discussed. 

The goal in the synthesis phase of .the research was not to 


recreate the.exact sound, but to.reasonably simulate it with .an.eye 


.to simplicity and minimum cost.. Practically, this goal was achieved 


with suitable timing circuits driving R-C oscillators.. As a check-on 


the validity of the analysis, this sound was also simulated by 


digital/analog methods. The waveform was mathematically described in 
the time domain, computed in discrete steps and converted to analog 
voltages. These voltages were then converted to sound energy. 

To show the practicality of synthesizing the bell, a device was 
designed and constructed which employed the bell sound as the warning 
sound for a vessel at anchor in reduced.visibility as required by 
U. 8. Coast!t! Guard Rules. To complete the suite of required fog 
signals for a sailboat, signals for sailing on various tacks and while 
under power were sagea,* This reduced.visibility warning device was 
designed for automatic signalling. Though the particular signals were 
for a sailboat, the concept is general enough for any small vessel 
where automatic fog signals may be desireable. Use of such a device 
on board small private, Naval, and Coast Guard vessels where the crew 


may be few in number and fully occupied with operating the vessel 


would be desireable. 


Refer to Appendix 6 for a discussion of the required warning signals 
for a sailboat in reduced visibility. 





SECTION 2 


THE CHARACTER OF A BELL TONE 


In the literature describing bell sounds or bell tones, the 
primary interest and discussions are related to the musical aspects of 
these natural sounds. Most good quality bells are described by their 
primary strike note in terms of the musical gcale. Past harmonic 
anlaysis of bells in acoustical research dealt almost exclusively with 
fine quality church or carillon bells. There doesn't appear to be 
much active research in this field today. During the period from 1920 
to 1935, considerable research effort was applied te.the problem. 

The physical explanation of the origin of sound from a bell is 
an extension of the notion of vibrational plates. The mathematical 
description of the flexural physics is beyond the intent of this 
chests, The exact mathematical solution of .this problem has not been 
obtained except for the case of thin walled bells. 

A bell after being struck gives off a sound composed of several 
separate frequencies (partial tones, or more commonly.— partials). 
Unlike many musical instruments which give off partials in a nearly 
harmonic ratio of 1:2:3:4..., a bell is not so constrained. An 
idealized ‘series for the bell partials would be 1.0#1.51:2.02:2.93: 
ee re An actual bell does not conform to-this.ideal series. 
The closer a bell is to this series, the purer is the note from a 


qualitative musical sense. 


* Refer to Lord Rayleigh's, The Theory of Sound, Vol I, p 388 
for a more complete treatise on the subject. 


** Authors disagree on this ‘idealized' eves lal] 
series given is 1.0:1.65:2.10:3.0:3.54:4.97... 


Another 
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There appears to be some disagreement among the various articles a 
written on the subject as to whether the strike note is generated by 
direct nodal vibrations, or if it arises as the result of a beat 


[4], [5] 


frequency. Curtis and Giannini appear to have employed a 
precise and controlled method for arriving at their results. They 
argue that the strike note from a particular bell which was analyzed in 
considerable detail arose as the result of harmonious blending of three 
close frequencies. 

The various partials of a bell have independent amplitude versus 
time response characteristics. And to make the wave form still more 
complex, the various partials may commence at separate times after the Py 
bell has been struck. The attack and decay rate of any single partial 
may be independent of all others. The higher frequencies of the 
composite wave form appear earliest in the spectrum after the bell is 
struck, and die away most rapidly. Some of the lower frequencies may 
not appear in the spectrum until as late as 1 to 3 seconds after the 
crash. The following table was extracted from a paper by Curtis and 
Stasaiai to illustrate the frequency content of a particular bell 
they studied. The bell chosen was a church bell tuned to the musical 
note F (345.3 hertz). 

The column headed 'Frequency of note' refers to the theoretical 
bell as a musician might describe it. The significant partials found 
to be present in that bell are in the column headed 'Tone Frequency’. * 
This illustration clearly shows that only the strike note came close 


to the desired frequency. Additionally, the strike note was straddled 


by other frequencies causing it to be in actuality a triplet. Curtis 
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® TABLE I 


Partials of a Bell Tuned to the Note F 





FREQUENCY TONE 
NUMBER PARTIAL OF NOTE NOTE FREQUENCY 
1 Hum note 172.6 Fi 160 
187 
2 Strike note 345.3 F 330 
345 
365 
2.4 Third (tierce) 410.6 Gi 385 
450 
3 Fifth (quint) 517.3 C 512 
4 Nominal 690.5 F 675 
e 700 
s 5 Upper Third 870.0 A — 
6 Wonder FeFeH 1034.6 ct 1060 





and Giannini went on to describe for this bell a total of some twenty 
significant partials and "a great many more partials which were weak 
compared to the ones which are recorded." 
To briefly summarize the complexities of the bell tone, one must 

list the following: 

1. The wave form is time limited. 

2. The frequencies that comprise the wave form may commence 

and end at independent times. 
3. The amplitude of each partial is generally time 
e independent of all others. 
4, The amplitude of each partial is time variant. 


5. The various frequencies present are generally not 


harmonically related. 
13 








SECTION 3 


RECORDING THE BELL 


The goal of recording the bell signature was to preserve a high 
quality data base from which to work. To achieve this end, it was 
desireable to have a good signal-to-noise ratio, preserve the relative 
magnitudes among the various frequencies present, and to obtain some 
flexibility in playback speed for sampling purposes. 

To minimize spurious background noise on the recording, the bell 
was placed in a large anechoic chamber and rung by an assistant. All 
recording equipment except for a microphone was placed outside the 
anechoic chamber. The microphone used was an Altec Lansing Model * 
21BR-150 broad-band microphone. The microphone was supported a distance , 
of 1 meter from the bell on a level with the bell's soundbow. 

The bell was suspended from its crown fitting and held rigidly in 
place so that no movement occured other than normal vibrational movement 
after being struck. The supporting structure for the bell was attached 
to ceiling and floor fixtures provided for this purpose in the chamber. 

Additional required equipment for the microphone outside the 
chamber was an Altec Lansing Model 526B microphone power supply. The 
response characteristics for the microphone and its associated power 
supply are typically within +1 db from 10 to 3000 hertz and within 
+3 db to 15000 hertz. A Hewlett Packard Model 466A broad band (DC to 


20000 hertz) amplifier with a selected gain of 20 db was used as a 


preamplifier prior to the tape recorder. Py 
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Block Diagram of the Recording Process 


The recording device chosen was an Ampex Model CP-100 instru- 
mentation tape recorder. This choice was made due to the excellent 
linear frequency response characteristics offered by the recorder 
over the anticipated frequency range of 100 to 15000 hertz. Two 
other salient features available on this tape recorder also contri-~ 


buted to its selection. These were, 


1. A frequency modulated recording and playback capability. 
This feature provides for a linear frequency translation 
without the necessity for amplitude compensation when 
the playback speed is different than the record speed. 

2. A wide choice of speeds for frequency translation 
purposes to add flexibility to the sampling procedure 
(this feature will be elaborated on in more length 
in the section on analysis). 


The signal was FM recorded at a speed of 60 inches per second. 
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This speed and recording method provided a tape recorder band-pass 4 
from essentially DC to 20000 hertz for the recorded signal. The overall 
band-pass of the recorded signal was 10 to 15000 hertz, limited by 
the microphone. 

Prior to recording the signal, the tape recorder was aligned so 
that the non-linearity did not exceed 0.75% (minimum achievable with 
the given equipment) over the recording-range.. An Ampex Model TC-10 
alignment set was used for aligning the tape recorder. The tape.heads 
on the recorder were cleaned and demagnitized prior to recording to 
insure a good signal-to-noise ratio.. Four different constant frequency 
test signals were placed on the tape recorder for calibration purposes - 
after analysis. The frequency of each of these signals was known to 
within 0.1 hertz. 

No attempt to measure the total power output of the bell was 
made due to the quite complex three dimensional. sound intensity pattern 
expected from the bell. This was-not required for the analysis, since 


the primary concern was to preserve the relative intensities among the 


various partials present in the bell signature. 





the Fourier Transform or Fourier Integral. It is a class of procedures 
for transforming a time series (discrete data samples) to its finite 
Fourier series. Many methods have been proposed and demonstrated over 
the years since Runge and Konig first described their duoeedunes, |e 
The history of modern techniques, expecially the Fast Fourier Transform 


[9] 


: SECTION 4 
ANALYSIS 
4.1. Discrete Method 
Discrete Fourier analysis simply stated is an extension of 

(FFT), are both interesting and well documented. The FFT, a special 
case of the discrete Fourier transform (DFT), is an algorithm for 

e seciesenly computing the finite series transform of the discrete data 
set. Its application is suitably fitted to discrete eomputation on 
digital computers. It finds wide application in digital spectral 
analysis, filter simulation, convolution.and other related fields. 
The significant feature that makes this clever. technique appealing over 
earlier techniques is the rapid method used to perform the desired 
operations. Time and money are interchangeable when discussing digital 
computation. For a comparison of computational time required for the 
FFT as opposed to earlier DFT methods, consider a time series consisting 
of N= 2" samples. To perform a discrete- Fourier transferm.using FFT, 
(N log, N) computational steps are required. For earlier more direct 
methods, ir computational steps would be required to perform the same 

‘ transformation . For a time series consisting of N = 1024, approximately 
a 50 to 1 savings in computational time is realized through use of the 


per, [10] 
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The discrete Fourier transform is defined by 
N-1 


A = xX exp (-zrirk/N) r = 0,1,...,N-1 (1) ‘ 
i 


4 
where A. is the rth coefficient of the DFT and Xx denotes the 
kth sample of the time series which consists of N samples and i = ¥-1 . 
The relationship between the DFT and the Fourier transform is shown 
in a paper by Cochran, Cooley, Favor and others. Since the FFT is an 
implementation of this definition, this relationship also defines 
the FFT. 

Given a time series with a constant sampling interval (@t) between 
each successive sample, the sampling frequency is given by 

f = 1/at we Baers (2) 


By the Nyquist sampling criterion, the resultant bandwidth recovered 


from the signal of interest would be * 
Band Width = B = 1/2At ) 
=f /2 (4) 
s 


This of course is only true if the signal is band-limited to B 
before being sampled. Given a frequency higher than B in the original 


(B < £, < 2B), then a finite Fourier series for the 


[12] 


signal, say f 


1 1 


signal after analysis would reveal a spectral line at £78. This 
aliasing can be observed with a stroboscope and a rotating machine, or 
more amusingly, by watching the wheels of a stage coach in a movie 


appear to rotate in the incorrect direction. To recover correctly 


the frequencies present below B, the signal to. be sampled must be 


filtered before sampling. Even with filtering, some error is introduced 





since it is impossible to completely band-limit a signal. Practically 
speaking this error is small and can be ignored if the filter chosen 
has a high roll-off and the corner frequency is chosen with care. 

The FFT yields a finite spectrum of N/2 distinct lines for N 


sampled data points. Therefore these lines will be separated by Af where 





» Boo. 228: 
- N/2 N 6) 
ts 
ae (6) 


With a time variant signal (the coefficients of the associated Fourier 
series vary with time) this simple relationship yields a paradox of 
accuracy. For a fixed sampling rate. (i.e., f is fixed), N must be 
made large with respect to f to recover a small Af. But if NAt 

is large compared to the time over which the coefficients vary signifi- 
cantly, then the coefficients recovered by analysis will be averaged 
over the time series duration. If N is made small with respect to f 
the average coefficients recovered will be closer to the true value at 
the beginning and the end of the time series. The penalty paid, of 
course, is that Af would now be larger and the spectral lines would 


[13] 


be farther apart. This paradox is yet another form of the well 
known uncertainty principle. A desireable compromise would be to have 
the coefficients change only a small amount over the period NAt: and 
yet have a sufficiently small Af to discriminate between adjacent. fre- 
quencies present in the spectrum. Stated in. another way, the goal in 
analyzing a non-periodic wave form is to achieve a quasi-stationary 


process over which Af and the resultant coefficients can be meaningful. 


There are many permutations of f and N to obtain such a result. One 
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technique used in this research was to analyze a time series once to 
obtain a fine Af, and then to analyze the same series again to obtain 
reasonably accurate coefficients. Ae 

For the spectral analysis of the bell tone, the recorded signal 
was filtered, amplified, sampled at constant intervals, and then stored 
on magnetic tape as the discrete.time series. The recorded signal. 
was played back for sampling at 1 7/8 inches per second which yielded a 
frequency translation of 32 to l.over the recorded speed.of 60 inches 
per second. This was done primarily due.to sampling rate limitations 
caused by the manner in which the sampling was performed. By sampling 
with the program as written (see Appendix 3), samples were written 
on magnetic tape after a set of 128 were collected. The upper sampling 
rate achievable by this method is limited. by. the. magnetic tape write 
speed, which is around 1800 hertz for the stated record length. The » 
sampling program stored 520 sequential records of 128 samples per record 
on a 7 track magnetic tape. The sampling. frequency used was 1024.0 
hertz. When this sampling frequency is translated by 32 (to correct 
for the tape speed reduction), a true sampling speed of 32,768 hertz 
is realized. By the Nyquist sampling criterion, this sampling frequency 
produces a bandwidth of 16,384 hertz. This band width is slightly in 
excess of that of the recorded signal. 

Prior to sampling, the signal. was amplified to a peak value of 
about 60 volts to minimize the noise introduced by the sampling process. 
Between amplification stages the bell tone was passed through a continuous 


band-pass filter. Ths band-pass of the filter was flat from 10 hertz to 


150 hertz and was down 3 db at 260 hertz. This upper 3 db point corresponds , 





to 8320 hertz when translated. This band-pass filter upper limit may 
seem low based upon the band-pass of other phases of the recording and 
sampling process. Later analog analysis showed that the original signal 
did not contain significant energy in frequencies higher than 6000 
hertz. 

The sampling frequency was chosen at 32,768 hertz (actual) since 
this is a power of two and corresponds to an integer separation of 


frequencies in the spectra for a sample size of N = 4096. By formula 5, 





Ag . 32768 
* ~ “4096 


8 hertz 

For complete analysis of one segment of the signal, the record 
size per window was chosen at 4096 samples. This corresponds to a 
record time length (NAL) of 125 ms. Each window overlapped the preceding 
window by 50 per cent. The time at which each window was analyzed was 
considered to be the time at the center of the window. Thus the coeffix 
cients from window 1 (time of window from 0.0 to 125 ms after the bell 
crash) were considered to exist discretely in time at 62.5 ms after the 
crash, and so forth for the remaining windows ef observation. A total 
of 32 windows were analyzed for this record size. 

A second complete analysis of the same time series was performed 
using a record containing 16384 samples per window. The time duration 
for this record length was 500 ms. The computed speetral line separation 
for these parameters is 2 hertz. Similarly to the first pass, a 50 per 
cent overlap of each successive window was employed. 

Sampling of the signal was performed on a hybrid analog-digital 


computer using an external frequency source for the sampling frequency 
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reference clock. The hybrid installation consisted of an SDS 930 
general purpose digital computer.interfaced with a Comcor CI-5000 
electronic analog computer. The analog-to=digital converter had an 
fourteen bit word length to represent discrete levels of 12 millivolts, 
based upon the analog variable range of +100 volts to -100 volts. 

The Fourier analysis computations were accomplished on an IBM 
360/67 digital computer (see Appendix 5 for the program used) since 
the program required for the record sizes employed exceeded the memory 
size of the SDS machine. Due to the different word sizes of the SDS 
and IBM machines, an assembly language subroutine was written to convert 
the sampled data word format. This subroutine is given in Appendix 4. 
The word format change was made on the IBM computer. For rapid selection 
of random windows from the entire time series, the sequential time 
series was stored on a pseudo-random access disk pack. Each pass of 
an analysis took the desired time series sub-set from the disk pack, 
analyzed the series, and then printed out a permanent record of the 
Fourier series coefficients. For selected portions of the analysis, 
graphs were drawn by peripherals to the IBM computer. 

The window size consisting of 4096 data points was considered 
sufficiently short to give reasonable aeeuracy to the resultant 
coefficients. The analysis using 16384 samples per window was used to 
determine the center frequency of the broadened spectral line for each 
coefficient. From these procedures, the frequencies stated in the 
following table are considered to be the significant ones in the 


original bell tone. 
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TABLE IT 


Partials Of the Fog Bell 


Significant Commence (End) Maximum Amplitude 
Frequencies (sec) (relative energy) 


565 0.125 1.6 0.24 
1370 0.125 1.5 0.15 
2331 0.0 0.875 1.76 
3061 0.0 2.0 1.05 
3320 0.0 1.0 1.68 
3f72 0.0 1.4 0.80 


The time listed when the partials commence are estimates since this 
information is relatively uncertain. The time when each partial ends 
is based upon the time when they fall to 0.01 per cent of their maximum. 
Many other frequencies were present in the spectrum, but these were 
either too short in duration or too insignificant in energy to analyze 
in detail. 

The frequencies 2331 and 3320 hertz contained the most significant 
amount of energy in the series for the bell. Since these frequencies 
commence early and die off fairly shortly after the bell crash, it is 
felt that these partials comprise the distinctive sound of the bell 
crash, It would appear that the frequencies 565, 1370, and 3061 give 
the bell its prolonged sound as it dies away. No attempt has been 
made to correlate the analyzed data with a musical scale or give a 
qualitative explaination of the bell sound. This. was not done since 
the bell chosen was for fog signaling and was.not tuned to-.any 
particular musical scale. The analyzed data conforms generally to the 


theory and format of the bells described by Curtis, Giannini and 


[14], (25) 


others. 





4.2 Continuous Method 
To check the results of the analysis by the discrete method, 

a continuous band-pass technique was employed. A special purpose 

audio spectrograph (Kay Missilyzer) was used for this task. This 

spectrograph records the signal to be analyzed on an endless magnetic 

tape which is mounted on a rotating drum. The spectrograph triggers 

a 5 ms integrator with a tuneable 20. hertz band-pass input on the same 

position of the drum each rotation. The magnitude of the output of 

this integrator is burned on a recording paper so the spectrum may be 

preserved. As the drum turns, the filter advances each rotation.of the 

drum. For the recording speed selected, the band width of the spectrum ° 

analyzed was 5000 hertz. The position on the drum where the integrator 

is triggered can be selectively pbsitioned so that successive slices ‘ 

(time windows) of the spectrum can be made. By manually transferring 

the spectral line amplitudes, a time. plot of amplitude versus time 


for the various partials can be developed. 


4.3 Comparison of Methods = 


The partials in the bell tone found by.the discrete method were 

also found to be present by the continuous method. The amplitude 

versus time information correlated between. the techniques fairly well. 

Since.the accuracy for the coefficients obtained by the continuous 

method should be greater due to the much shorter window size, one 

would not expect the amplitude versus time plots to match exactly.. = ’ 
By using 16384 discrete data points, the accuracy of the frequencies 


of each partial is + 2 hertz. No theory is known to the author to 
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develop the bounds on the accuracy of the coefficient magnitudes for 

the continuous method. It is considered.that this amplitude information 
is quite dependable. Appendix .2 shows the amplitude versus time. plot 
for each of the six significant partials found through the discrete 

and continuous methods. 

Using a lower sampling frequency for the discrete analysis 
method, the bandwidth of the time series could have been reduced. If 
this were done, it would have been possible to use a shorter window 
length and still maintained a small Af. The result of the analysis 
would have been amplitude information with higher confidence. 

One distinct advantage of the discrete method over the continuous 
one for some applications is an ability to present phase information 
about each frequency present in the.spectrum. This additional piece 


of information was not required for this research since the ear cannot 


[16] 


determine phase information about a complex wave form. 





SECTION 5 


REDUCED VISIBILITY WARNING EQUIPMENT 


5.1 Discussion 

Quite often aboard a private sailing yacht, a crew may be 
sufficiently occupied with sailing the craft during reduced visibility 
that sounding fog signals could. be.a burden. Sounding such repetitive 
timed signals is a boring but. quite important. task. Additionally, the 
specified ‘at anchor' warning signal may be required at a time when no 
crew members are on board the craft. A semi-automatic signalling device 
could alleviate the problems created by the aforementioned examples by 
providing reduced visibility warning signals for the craft. Such a device 
ideally should be simple to operate while providing dependable contin- 
uous warning signals for underway and at anchor operations. 

The next two sub-sections describe the design for a device 
which can generate these signals. A later section describes the 
synthesis of the bell tone which was included in the device. The final 
form of the device was constructed out of solid state devices and placed 
on printed circuit boards. This device could be packaged as a small 


portable unit or permanently mounted aboard a yacht. 


5.2 Timing Circuitry 

The low visibility warning. equipment.has basic periodic 
features. These are dictated by the methods of generating the various 
sounds and the requirements for these sounds. . 

The major periodic feature common to the anchored warning signal 


and all the low visibility warning sounds while underway is their period 
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between sounds. . All of these sounds (with the exception of the required 
signal for underway under power in international waters) have a maximum & 
period of one minute. The maximum specified period for under power in 
international waters is two minutes. Prudent seamanship dictates that 
these intervals not be fixed. over any long. period of time. This 
desired aperiodicity prevents two.vessels from sounding similar warning 
signals synchronously. Synchronism.could bring about the. undesireable 
side effect of a collision and thus the. ruination of one's day. 

A master clock is required for the warning equipment and is 
specified by a frequency of from.less than one cycle per minute to 
less than one per two minutes... Unijunction oscillators are immediately 
suggested by their extreme simplicity.and ability to satisfy these * 
requirements. The basic form of .a-relaxation oscillator with a periodic 
pulse output was chosen and is shown.in Figure 1. A potentiometer is 
used to vary the oscillator frequency manually when desired. With the 
component values shown, the periodic range is variable from approximately 
55 seconds to 90 seconds. 

The periodic pulse output of this master oscillator was fed to. 
the trigger input of the first .of a series.of five mono-stable (one-shot) 
multivibrators. These one-shots serially generate all.the required 
logic levels.for sounding the various warning signals. The pulse from 
the master oscillator turns on.the first one-shot. When the first one 
..turns off, the second one-.turns on, etc. This turn on/turn off procedes 
through all five one-shots until the last.one is off. This sequence ‘ 
is initialized by each trigger pulse from_the master oscillator. . 


The turn on of the next one-shot is accomplished by inverting the previous 


28 








O3.6 V 
th 
- 

6 
ONE-SHOT FROM 


yL9l4 AND DISCRETE 
COMPONENTS R3,CI, AND C2 


FIGURE 3 






TRIGGER 
INPUT 


pulse and then differentiating it. Direct differentiation of the 
‘conjugate’ wave-form of the previous one-shot is impractical since the 
trigger pulse propagates through on.this wave-form. Additionally, 

some isolation is required for the devices chosen since drive capability 
is limited. 

The 'on' time pulse of each one-shot is shown in Figure 2. 
Selection of the proper pulse train for the various sounds is accom- 
plished by a manual function selector switch. This switch is in 
actuality a variable input AND gate. 

The active devices chosen for the one-shots were Fairchild uL914 
dual two-input NAND gates. By the use of one external resistor (R3), 


[17] Figure 3 shows 


and capacitor (Cl), these gates become one-shots. 
the internal circuitry of these gates and the application of the external 
components. These devices are designed.for high-speed digital logic 
-applications and as such turn.on.with. small signal levels. Spurious 
noise and small supply voltage variations can cause unwanted triggering 
of the one-shots. To minimize this occurence, a large. capacitor, C2, 
(typically, 50 uf or greater) was placed.as shown in Figure 3. 

These micrologic devices employ a supply voltage in the range 
of 3.0 to 4.2 volts. The nominal recommended voltage is 3.6 volts. 
A zener diode was employed to give this desired supply level. A 1000 
wi capacitor was required in parallel with.the zener to give additional 
stability to. the supply-voltage. Without this capacitor, spurious 
triggering results and the chain of one-shots frm an oscillator. With 
the capacitor, some spurious triggering still results, but the one-shot 


chain does not go into continuous oscillation. This spurious triggering 
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could also be eliminated by providing a constant-voltage, variable- 
current power source. This further complexity and cost is not required 
in this application since the correct output from the timing chain 

is achieved. 

After the function switch selects the proper chain of pulses, 
these pulses are applied to a relay actuated switch which connects the 
supply voltage to either the horn oscillators or to the bell circuitry. 
The supply voltage for the bell oscillators also goes through a 
transistor switch which forms the voltage wave form shown in Figure 4. 

The pulses from the chain.of one-shots is also applied to a NAND 
gate. This gate develops the logical voltage that switches in a 
listening section and turns off the power amplifier when signals are 
not being sounded. This technique of listening is similar to a simple 
intercom system. The output speaker for the power amplifiet acts 
as the microphone input to an amplifier-speaker combination during this 
listening period. This remote listening device provides a degree of 
safety for the passengers and crew of ones own vessel during periods of 
reduced visibility. This listen feature can also be selected continu- 


ously by the master function selector switch. 


5.3 Description Of The Entire System 


Figure 5 shows in block diagram form the relationships of the 
various sub-parts of the entire reduced visibility warning equipment. 
Functionally, the system provides.the following features: 


1. The three required reduced visibility signals for a 


sailboat underway. 








FIGURE 5 REDUCED VISIBILITY WARNING 


EQUIPMENT 
-SIMPLIFIED BLOCK DIAGRAM— 







OUTPUT 
SPEAKER 





LOUD-HAILER 


LISTEN 





2. A reduced visibility signal for a vessel underway under 
power. » 
3. A manually operated warning horn for entering or leaving 
blind channels and slips. 
4, An automatic and manual listening device. 
5. A loud hailer. 
6. A rapidly ringing bell for a vessel at anchor in reduced 
visibility. 
7. Test positions for 2 and 6 so these functions may be 
checked out in port at reduced volume for preventative 
maintenance. 
The final power amplifier common to all features except 4 is a * 
standard class B transistor power amplifier. This output stage, the 
oscillator for all underway signals, and the loud hailer circuitry were 
taken directly from a commercial fog horn/loudhailer device. 
The underway signals employ a unijunction oscillator of a nominal 
frequency of 200 hertz. This is either actuated by the selected timing 
chain associated with a specific signal (1 or 2) or manually activated 
by a push-button switch for feature 3. When activated manually, a 
different load resistor is placed in the oscillator which causes 
an output frequency of nominally 380 hertz. 
The listening device partially described under the timing 
circuitry section employs a commercially available 500 milliwatt 
direct-output audio amplifier and an 8 ohm water proof speaker. The 7 
output speaker for the power amplifier is switched to act as a micro- 


phone input to the listen amplifier. A 3-pole double-throw relay 


34 


| eee 


actuated switch is used to remove power from the output stage to 
prevent damage to the transistors, and to switch the speaker to the 
listen amplifier input. The manual selection of 'listen' is achieved 
by placing a constant drive voltage on.the base of the relay driver. 
The automatic ' listen' feature is accomplished by placing the output 
of a NAND gate as described under.the timing circuitry section on the 
base of the same transistor. When installed on. a vessel, the ‘listen' 
output and the calling and emergency channel of a ship-to-shore receiver 
could be mixed and placed on a speaker.in the cockpit near the helmsman. 
This entire system provides fer a.typical sailboat all the 
advantages of semi-automatic warning.and signalling devices. The 
object of the equipment is to provide greater safety and flexibility 


for yacht sailing, motoring, or anchoring in a reduced visibility 


environment. 





SECTION 6 


SYNTHESIS OF THE BELL TONE 


cree Synthesis By Discrete Computation 


The bell sound was described mathematically as the superposition 
of the six most prominent partials found during the analysis. The 
partials were written as time variant sinusoids. The time response of 
each partial was approximated by fitted exponential curves. 

A FORTRAN II language program was written for an SDS 930 digital 
computer for computational purposes. The equivalent of sampled data 
was computed in discrete intervals of .000125 seconds. The computational 
step size in seconds and the frequency present in the wave form were 
based upon a 'sampling' frequency of 8000 hertz. 

A total of 8000 samples were computed and stored in a data table 
for a record length of 1 second at the stated clock frequency of 8000 
hertz. Provisions were made in the program to permit parameters of 
attack, decay, and amplitude to be varied. These provisions were made 
so that some experimental modifications to the wave form could be made 
in light of qualitative analysis. 

A machine language (META-SYMBOL) subroutine callable by FORTRAN II 
was written to perform the task of digital-to-analog conversion. This 
subroutine was controlled at a rate determined by a clock on the 
associated analog computer. If the clock were something different from 
the program oriented 8000 hertz, the equivalent of frequency/time 
translation would be performed on the data. The resultant analog 
voltage from this program was passed through a pass-band filter with a 


band width from 20 to 4000 hertz to minimize sampling noise. The wave 
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form was then amplified and reproduced through a speaker. The basic 
program had provisions for continuously repeating the same one second 
data record. The resultant effect was a bell being struck at one 
second intervals and ringing until being struck again’ 

The results of this experiment were quite encouraging. The 
bell-like sound that resulted was considered to be a reasonable likeness 
to the original sound. Therefore, the decision was made to proceed on 
the assumption that simple exponential approximations to the time 
response of each partial would be satisfactory for an engineering 
approximation. It was felt that some experimentation with the rise and 
decay times and the maximum magnitude of each partial would be necessary 
to optimize the sound. 

It was felt that this synthesis technique was a useful tool in the 
overall project as a verification of the engineering assumptions for the 
synthesis and analysis techniques used. If an investigation of more 


complicated sounds such as those involving voice inflection and accents 


were being conducted, then this step would have been invaluable. 


6.2 Synthesis of the Bell Tone by Electronic Circuitry 


Any analyzed sound can be exactly duplicated by man if sufficient 
complexity of circuitry and design time are expended. Practically 
speaking, this exact duplication is seldom desireable. Basic engineering 
concepts dictate that some of the objectives to be pursued when designing 
a portable sound production device are that it should be small, light-— 


weight, relatively inexpensive, reliable, and require a small amount of 


power. Looking at the desireable aspects of a practical bell simulator, 





one would expect it to generate a reasonable bell-like tone from a 
small inexpensive device. Hopefully not all the partials present in 
the original bell sound would be required for reasonable simulation, 
and the required partials could be simply generated. 

The first technique considered for synthesis of the bell tone 
was the use of a very high Q band-pass circuit (400 < Q < 1000). Such 
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circuits are possible using active filter techniques. Such a 
filter would employ only active devices, capacitors, and resistors. By 
causing such a filter to ring by introducing the required driving 
function, it would be possible to have a rising and decaying sine wave. 
By superposition of several of these waves, it would be possible to 
build up the synthetic bell tone. To achieve these very high Q 
circuits, a multiple-pole filter is required which would dictate 
several active devices per filter. As a further disadvantage, the 
rise and decay rates of the sine waves generated in this manner would 
not be independent of each other. 

The second technique investigated was the use of a constant 
amplitude sine wave oscillator fed into a gain-controlled amplifier. 
By varying the gain of this amplifier, the resultant wave shape would 
simulate one partial. Some wave shaping circuitry would be required to 
control the gain of this amplifier. There are several inexpensive 
integrated circuits presently being marketed that could be employed as 
the gain-controlled amplifier. These integrated circuits are generally 
designed for intermediate frequency (IF) amplification and they have 


automatic gain control (AGC) circuitry built into them. Typically, an 


80 db control range can be achieved for a few volts of AGC voltage. 





The manufacturers of these devices list the useful frequency range of 
these devices from DC to several megahertz. This method of tone 
generation was put aside when a still simpler method was discovered 
and investigated further. 

The method finally chosen to generate the individual partials 
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was twin-t oscillators driven by a variable supply voltage. 
significant feature of the twin-t oscillator for this application is 
that it can be made sufficiently frequency stable during variations 
of supply voltage. Most oscillators vary in output frequency and 
amplitude as the supply voltage varies. For the twin-t oscillator, 
the frequency variations can be kept to less than 1 1/2 per cent at 
the design frequency for large (0 to 10 volts) excursions of supply 
voltage. The resultant output magnitude varies almost linearly with 
the supply voltage. 

The oscillator is made frequency stable by selecting a high- 
current gain transistor (he, from 150 to 200) as the active device 
and including a large resistance in the base feedback path. If the 
supply voltage has a wave form as shown in Figure 6, then RC form 
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a load charging time constant during the time T.: During the time Typ 


CRA form a discharging time constant. The resultant supply voltage 


as seen by the oscillator is given by the following equations. 


v(t) = Vb (1 - exp(-at)) O<t <1, (7) 
On = R. Ce (8) 
v(t) = V(T,) exp(-a,t) T<£tat (9) 
o = R,C (10) 


d d f£ 
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It is implicitly assumed that the partial amplitudes versus time 
can be so simply approximated. As it turns out, this assumption is 
correct, the approximations of the partial amplitudes versus time being 
non critical. The true test is not the mathematical justification, but 
rather the ear of the observer. Compromise must be made for the goal 
of simplicity and low cost. 

The choices of Rus Ry and Ce affect both the amplitude (R. and 
Ry comprise the load) and to a lesser degree the frequency of the 
output. Variations of Ry will cause the oscillator to tune over nearly 
an octave. For optimum frequency stability, the oscillator should be 
adjusted to the center region of its tuneable range. Convenient design 
thumb rules are C, = 2C, , and R, = 0.1 Ri3 the values for C, and R 


2 1 2 1 1 
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can be selected from easily used nomograms. The resistor Ro is 
used as a feed to a mixing bus with other oscillators. Ro was also 
used to limit the gain of each oscillator to the desired value and 
provide a degree of isolation from other oscillators. The mixing bus 
was the input to an emitter follower amplifier which provided a high 
input impedence for the oscillators. 

For the first attempt at synthesis, six oscillators corresponding 
to the six most prominent partials (those at 565, 1370, 2331, 3060, 
3320, and 3773 hertz) found in the original bell tone were constructed. 
Three of these frequencies (2331, 3320, and 3773 hertz) reach their 
peak amplitudes earliest and probably give the bell its distinctive 


sound shortly after being struck. An attack time for all of these 


oscillators was chosen at 100 milliseconds (ms). The remaining three 


frequencies reach their peak amplitudes significantly later and 








probably give the bell tone its hum note. These latter frequencies 
were grouped together and given an attack time of 300 milliseconds. 
The values chosen for attack times for each group corresponds roughly 
to the average attack time for that group found in analysis. Admittedly, 
this broad grouping is an oversimplication, but certainly desireable 
from the view point of minimizing the required circuitry. 

The supply voltage for these oscillators was generated by employing 
a grounded emitter astable multivibrator driving transistor switches 
(see Figure 4). The on time of these switches corresponds to the attack 
time of the partials. The periodic nature of the astable multivibrator 
yields the effect of a bell being repeatedly rung. The wave form of the 
astable multivibrator was a rectangular wave having an on time of 100 ms 
and a period of 700 ms. For the longer attack-time group of partials, 
a monostable multivibrator was used to extend the on time of the 
associated switch to 300 ms. A relay actuated switch in series with the 
transistor switches was used to control the on time of the ringing bell 
to 5 seconds out of every minute. 

By properly adjusting the maximum amplitude of each oscillator, 
the resulting bell tone was a fair representation of the original sound. 
The rapid periodic ringing of the bell suggested that some of the lower 
frequencies comprising the longer lasting hum note group could be 
eliminated. Hopefully this could be done with little degradation of the 
simulated bell tone. By experimentation, it was found that only the 
partials at 2331 and 3320 hertz were required to represent a bell being 


repeatedly rung at intervals of 700 ms. The final form of the bell tone 
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simulator therefore contained only two variable amplitude oscillators. 
In light of the stated objectives, it is felt that the simulated bell 


tone adequately fulfills the requirements for a warning device of an 


anchored vessel. 





SECTION 7 


SUMMARY 


7.1 Analysis 


Both the discrete and continuous analysis techniques suffer certain 
limitations. These limitations arise due to the finite frequency 
spectrum representation of time and band limited complex wave forms. 

The resulting inaccuracies are inevitable and arise as a result of the 
uncertain relationship of the time signal to its frequency transform. 

By selection of proper sampling frequency and time series duration, 
a quasi-stationary process can be approached for the purpose of obtaining 
both frequency and amplitude information. Permutations of sampling 
frequency and time duration of the series allow optimization of frequency 


information, or amplitude information, but not both. e 


7.2 Synthesis 


Complex wave forms can be approximated quite simply using solid 
state devices if the exact wave form is not to be duplicated and 
compromises can be accepted. This synthesis can also be done in the form 
of a mathematical model if desired, so the model can be examined more 
closely before expensive and time consuming 'bread-boarding' is attempted. 

It is felt that the low visibility warning equipment has demon- 
strated the feasibility and practicality of such a device. The object 


was not to develop a revolutionary apparatus, but rather to demonstrate 


the engineering techniques and expertise for such a device. 
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APPENDIX 1 


BELL SPECTRUM BY DISCRETE ANALYSIS 


Appendices 1.1 through 1.6 show contiguous time windows which 
have been transformed to the frequency domain for the first 750 ms of 
the bell tone. Of interest is the rapid fall of the spectral line 
amplitudes above 4000 hertz, and the slow rise of the lower frequencies. 


For these spectra, a sample size of 4096 data points was used for the 


transformation which yielded a Af of 8 hertz. 
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APPENDIX 2 . 


COEFFICIENT AMPLITUDES VERSUS TIME FOR THE BELL 


By discrete and continuous analysis, the magnitudes of the 
coefficients for each significant partial was found as discussed in 
Section 4. Appendices 2.1 through 2.6 show a comparison of these 
coefficients as a function of time. As discussed in Section 4.3, the 
plots would not be precisely the same. 

For those coefficients with a slow rate of change (e.g., 565, 
1370, and 3061 hertz), the results of the two methods used compare 
fairly well. For the partials with a faster time rate of change 


(particularly 2331 and 3320 hertz), the correlation of the results 


was poorer. The general shape of these plots were similar however. + 
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APPENDIX 5 


FAST FOURIER TRANSFORM ANALYSIS PROGRAM 
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APPENDIX 6 


REQUIRED FOG SIGNALS FOR A SAILBOAT 


Definitions for the sailing terms used: 


Relative wind 
in this sector 
is a starboard tack 


Relative wind 
in this sector 
is a port tack 








Relative wind in this sector 
is called wind abaft the beam. 


Sound signals for sailboats in reduced visibility while 
underway in Inland and International waters: 


Operation of Signal Maximum interval 
sailing vessel between signals 
(minutes) 
Starboard Tack 1 blast of 1 
foghorn 
Port Tack 2 successive blasts 1 
of foghorn 
Wind abaft the beam 3 successive blasts 1 
of foghorn 
we* 
While motoring 1 prolonged blast 1 (Inland) 
of foghorn 2 (International) 


A sailboat shall sound at intervals not to exceed 1 minute 
a rapidly ringing bell for about five seconds. 


* Blast is defined as a duration of not over 2 seconds. 


** Prolonged blast is defined as a duration from 4 to 6 seconds. 
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APPENDIX 7 
“ SCHEMATICS FOR THE REDUCED VISIBILITY WARNING 
EQUIPMENT 
This appendix contains an expanded block diagram for the reduced 


visibility warning equipment and schematics for the various sub-parts 


which the author designed. 
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